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Abstract 



Different phases of sohd boron under high pressure are studied by first 
principles calculations. The a-Bi2 structure is found to be stable up to 270 
GPa. Its semiconductor band gap (1.72 eV) decreases continuously to zero 
around 160 GPa, where the material transforms to a weak metal. The metal- 
licity, as measured by the density of states at the Fermi level, enhances as 
the pressure is further increased. The pressure-induced metallization can be 
attributed to the enhanced boron-boron interactions that cause bands over- 
lap. These results are consist with the recently observed metallization and the 
associated superconductivity of bulk boron under high pressure (M.I.Eremets 

et al, Science293, 272(2001)). 
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It is well known that many semiconductors and insulators undergo nonmetal-metal transi- 
tions upon applying high pressure. Ehi Most of these metallization are associated with a struc- 
tural transition from low coordination insulator to a high coordination metallic phase.Bi^ In 
some case, the pressure-induced metallization can be attributed to the increased interatomic 
interactions which lead to the band overlap.il Recently, pressure-induced metallization and 
superconductivity were found in solid boron above 160 GPa.i There is no direct experimen- 
tal evidence whether the nonmetal-metal transition is caused by the band gap closure or a 
structural transition into the metallic bet phase.0 Moreover, the superconducting transition 
temperature increases from 6 K at 175 GPa to 11.2 K at 250 GPa. This is in contrast 
to typical sp metals such as Al, where the Tc decreases with pressure. i'i In this letter we 
present results of first principles calculations on solid boron under high pressure for different 
phases (a-Bi2, bet, and fee). Our results suggest that a-Bi2 phase is stable and maintains 
its rhombohedral structure up to 270 GPa. Around 160 GPa, a nonmetal-metal transition 
is found in a-Bi2 phase. Further increases in pressure leads to an enhanced density of states 
at the Fermi level. Charge density analysis shows that the pressure- induced metallization 
is mainly due to the increased interatomic interactions that broaden the band width. We 
propose that the recently observed metallization! can be understood in terms of the band 
gap closure. 

Under ambient pressure, solid boron exists in various complex crystal structures with 
icosahedral B12 cluster as a common structural component. The B12 icosahedrons can be in- 
terlinked by strong covalent bonds in different ways to form different polymorphs such as the 
a-rhombohedral B12 (a-Bi2), a-tetragonal B50, or the /?-rhombohedral B105 (/?-Bio5).0'0lli 
Among these crystal structures a-Bi2 has the simplest form with one B12 cluster per rhom- 
bohedral unit cell. Previous experimental studies on solid boron include the structural 
analyses^, bulk moduluJiil, Raman spectrum^, and electron density distributionlii. On the 
theoretical side, the structural, electronic and vibrational properties of solid boron (mostly in 
a-Bi2 phase) have been calculated by first principle methods.00'0Ei In particular, a struc- 
tural transition from the insulating a-rhombohedral to metallic body-centered-tetragonal 



(bet) structures was predicted to occur at the pressure of 210 GPa.EI But no detailed analy- 
sis on nonmetal-metal transition were done. 
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FIGURES 




FIG. 1. Optimized a-Bi2 rhombohedral crystalline structure shows strong bonding between 
icosahedral B12 clusters. Present GGA calculations give the cell constant a=4.98 A and cell angle 
a=58.2°, which compare well with experimental values of 5.06 A , 58.2°0. (See Table I for the 
detailed structural parameters.) 

Our density functional calculations is based on the plane-wave pseudopotential 
methodBS. The density functional is treated by the generalized gradient approximation 
(GGA)lli with the exchange-correlation potential parameterized by Wang and PerdewEl. We 
use a ultrasoft nonlocal pseudopotential to model the ion-electron interactions^. The kinetic 
energy cutoff for the plane-wave basis is 240 eV. Large set of Monkhorst-Pack k meshes0 ( 
8x8x8, 14x14x14, and 14x14x14 for a-B^, bet, and fee lattice respectively) are used to 
sample the Brillouin zone. For a given external pressure, both the cell parameters and the 
atomic positions are fully optimized. 

The optimized crystalline structure of a-Bi2 is shown in Fig.l. It is a network of B12 clus- 
ters with two kinds of intercluster covalent bonds {di and ^2 in Table I) between neighboring 
icosahedrons. The calculated cohesive energy 6.95 eV/atom is comparable to the experi- 
mental value 5.81 eV.ii The indirect gap 1.72 eV from our GGA calculation agrees with the 
estimated experimental value (2.0 eV)ii better than the previous IDA resuh (1.43 eV).0 
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Table I summarizes the structural parameters of a-Bi2 under different pressures. Experi- 
mental and previous theoretical results are included for comparison. Satisfactory agreement 
is found between present calculations and previous works. As shown in the table, upon 
relaxation the a-Bi2 maintains its rhombohedral lattice structure up to high pressure. The 
cell angle a between the lattice vectors does not change with pressure, while the lattice 
constant smoothly decreases. All the interatomic bonds are compressed with most notable 
changes in the long intercluster bonds (^2)- This leads to a dramatic modification on the 
electronic states (to be discussed later). 
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TABLES 



Table I. Summary of structural parameters of a-Bi2 under different pressures (0 ~ 400 
GPa, given in bracket), a denotes the lattice parameter; a describes the angle between the 
rhombohedral lattice vectors; r is the mean radius of the B12 icosahedron; di and 62 are the 
short and long intercluster bond lengths. Available experimental valueJiSlini and previous 
theoretical resultslii'0 are included for comparison. 





a 


a (A) 


r (A) 


di (A) 


d2 (A) 


Exp.ll2l (0 GPa) 


58.2° 


5.06 


1.69 






Exp.§ (0.5 GPa) 


58.0° 


5.07 


1.70 


1.67 


2.00 


Theor.0 (0 GPa) 


58.2° 


4.98 


1.68 


1.65 


1.98 


Theor.0 (0 GPa) 


58.1° 


5.03 


1.69 


1.67 


2.00 


This work (0 GPa) 


58.2° 


4.98 


1.70 


1.67 


1.99 


This work (100 GPa) 


58.0° 


4.56 


1.56 


1.51 


1.74 


This work (200 GPa) 


58.1° 


4.34 


1.49 


1.43 


1.63 


This work (300 GPa) 


58.2° 


4.19 


1.44 


1.38 


1.56 


This work (400 GPa) 


58.3° 


4.06 


1.41 


1.34 


1.51 
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FIG. 2. Relative enthalpy of bet and fee phase with respeet to a-Bi2. The crossing of bet and 
a-Bi2 phase oeeurs at about 270 GPa. No cross between fee and bet phases is found up to 400 
GPa. Insert: cohesive energy per atom of Q!-Bi2, bet, fee phases as a function of atomic volume. 

To investigate the pressure-induced structural transition we optimized both the cell pa- 
rameters and atomic positions for q;-Bi2, bet, fee phases under each external hydrostatic 
pressure. Both the total energy and the enthalpy for each phase were calculated up to 400 
GPa. The results are shown in Fig.2. At the ambient pressure, we find 0.82 eV energy 
difi^erence between the equilibrium q;-Bi2 and the bet phase, and 1.45 eV difi^erence between 
the q;-Bi2 and the fee phase. Upon applying pressure, q;-Bi2 remains the low energy struc- 
ture but the enthalpy of bet phase becomes lower than that of q;-Bi2 at 270 GPa. Such a 
crossing indicates a thermodynamic instability of q;-Bi2 and a possible structural transition 
to the bet phase. The enthalpy of the fee phase crosses that of the q;-Bi2 at the pressure of 
315 GPa, but it is never lower than that of the bet phase up to the highest pressure (400 
GPa) studied. 

Previous LDA calculations predicted a structural transition from q;-Bi2 to bet at 210 
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GPa.Q We have also preformed independent LDA calculations and obtained a similar tran- 
sition pressure of 200 GPa. Our GGA transition pressure is about 60 GPa higher than the 
LDA results. As GGA calculations provide a better agreement with experiments in the band 
gap than LDA calculations, we believe that the GGA transition pressure (270 GPa) is more 
reliable. In the previous LDA calculations, a further transition from bet to fee was predicted 
to occur at 360 GPa.0 We have not observed such a transition up to 400 GPa in both LDA 
and GGA calculations. The discrepancy may due in part to the fact that our calculations 
are fully optimized in both the lattice parameters and the atomic positions, while the c/a 
ratio of the bet structure was fixed in the previous calculations.0 



FIG. 3. The electronic density of states of a-Bi2 under different pressures. Gaussian broadening 
of 0.05 eV is used. The broadening of both the valence and the conduction bands with pressure 
leads to the band gap closure and the formation of metallic states around 160 GPa. 

From above discussions we conclude that it is unlikely that the experimentally observed 
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nonmetal to metal transition at 160 GPaB is due to the structural transition into the metallic 
bet phase. The alternative explanation for the metallization is the band gap closure upon 
external pressure, similar to that happened in iodine-S As shown in Table I, upon applying 
pressure, the a-Bu maintains its rhombohedral lattice structure with compressed intercluster 
and intracluster bonds. The strengthened boron-boron interactions broaden both the valence 
and conduction bands, eventually leads to the band gap closure and the development of 
metallic states near the Fermi level. This effect is clearly shown in Fig. 3, where we present the 
density of states for several different pressures. We perform further analysis of the electron 
density. Upon applying pressure, both the two-center intercluster bonds and the three-center 
intracluster bond&0 are greatly strengthened while the two-center intracluster bonds are 
less sensitive. This can be seen in the contour plot of the total charge distributions for 
different pressures. Fig.4 shows an example of the contour plot, which includes the two- 
center intercluster and intracluster bonds. A much higher charge density in the area of the 
intercluster bonds is found at higher pressures. Similar effect was also observed in the three- 
center intracluster bonds. Thus, we conclude that the enhancement of both intercluster and 
three-center intracluster bonding leads to pressure-induced metallization. 
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FIG. 4. The total charge density contour plot for a-Bi2 phase under different pressures: GPa 
(upper), 160 GPa (middle), 240 GPa (down). The current slice includes two-center intercluster and 
intracluster bonds. The representative sites for boron atoms are labeled as B in upper plot. Red, 
yellow, green, blue colors indicate decreasing electron density. A significantly enhanced charge 
density (red color) in the area of the intercluster bond can be seen in the case of higher pressures. 

A clearer picture of nonmetal-metal transition can be seen in the pressure dependent 
of the band gap shown in Fig.5. From extrapolation, the band gap closure occurs around 
160 GPa. Experimentally, it was found that /3-Bio5 undergo a nonmetal to superconductor 
transition at about 160 GPJ. Although our calculations were based on a-Bi2 phase, we 
believe that our results are qualitatively comparable to the experiments because both phases 
have the same icosahedral structural component. 
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FIG. 5. The band gap (filled dots) and the electronic density of states N{Ef) (open squares) at 
the Fermi level for a-Bi2 as functions of pressure. The extrapolated pressure for band gap closure 
is 160 GPa. In the metallic phase, N{Ef) increases rapidly with pressure at first, then tampers 
off around 240 GPa. 

In the metallic phase of q;-Bi2 our calculations show that density of states N{Ef) at the 
Fermi level increases rapidly with pressure from 160 GPa to 240 GPa, beyond which it sat- 
urates (see Fig.5). The enhancement of metallicity with pressure can be directly associated 
with the increase of superconducting transition temperature Tc from 6 K at 175 GPa to 11.2 
K at 250 GPa.i In the phonon mediated superconductivity, the Tc can be obtained as@'@ 

Tc = ri^"^(-^-0 ^-^'-0^62A;z- )' 

where is the characteristic phonon frequency, A is the electron-phonon coupling strength, 
fi* is the screened coulomb repulsive interaction. The electron-phonon coupling A is related 
to the density of states at the Fermi level N{Ef) by 

A = N{Ef) < > /{M < 00^ >), 

where < > is the average square of the electron-phonon matrix element and < cu^ > 
is averaged phonon frequency.S Typically for sp metal such as aluminum, < cu^ > in- 
crease with pressure and N{Ef) decreases, while fx* is insensitive.^ Thus, Tc decreases with 
pressure.!'! On the contrary, in solid boron, Tc was found to increases substantially with 
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pressure (175~250 GPa).Q Our calculations show that such increase in Tc can be understood 
by the dramatic increase in the density of states N{Ef) with pressure after the metallization 
of a-Bi2. In addition, we also have calculated the N{Ef) and its pressure dependence for the 
metallic bet phase. The N{Ep) for bet phase is much higher than that of a-Bi2 phase and it 
does not show a clear pressure dependence. This result further supports our interpretation 
that the observed metallization in boron may not originate from a structural transition from 
a-Bi2 to bet phase. 

In summary, we have performed first principles calculations to investigate the structural 
and electronic properties of solid boron under high pressure. Structural transition from a- 
Bi2 to bet phase is found to occur at about 270 GPa. The semiconductor gap (1.72 eV at 
ambient pressure) decreases upon external pressure. Band gap closure and a nonmetal-metal 
transition is found around 160 GPa. As the external pressure further increase (160 GPa to 
240 GPa), the density of states at Fermi level increase dramatically. We conclude that the 
experimentally observed metallization in boron might be due to the band gap closure instead 
of a structural transition into the metallic bet phase. 

This work is supported by the U.S. Army Research Office (Grant DAAG55-98- 1-0298). 
We acknowledge computational supports from the North Carolina Supercomputer Center. 
*: zhaoj@physics.unc.edu 
^: jpl@physics.unc.edu 
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